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1.  Executive  Summary 

The  DARPA  Topological  Data  Analysis  program  has  supported  the  development  of  topo¬ 
logical  methods  for  analyzing  large  sets  of  points  in  metric  spaces  [1-5].  The  goal  of  this 
project  is  to  develop  similar  methods  for  data  sets  on  which  the  natural  structure  is  not  a 
metric,  but  rather  a  partial  order.  These  methods  are  being  developed  for  application  to 
Boolean  data,  i.e.,  sets  of  bit  strings,  and  to  geometrical  data,  e.g.,  points  in  Lorentzian 
manifolds. 

2.  Research  Summary 

During  the  first  quarter  of  Year  1  of  the  grant  I  began  formulating  a  framework  for  topo¬ 
logical  analysis  of  points  in  Lorentzian  manifolds,  partially  ordered  by  causality.  This  for¬ 
mulation  was  the  basis  for  two  talks  I  gave  at  the  European  Marie  Curie  Research  Training 
Network  on  Random  Geometry  conference  Causal  Sets  2006 :  “Dimension  of  causal  sets' 
and  “Algebra  and  topology  of  partial  orders'1. 

During  the  second  quarter  I  developed  definitions  for  filtered  topologies  on  antichains  in 
partial  orders  based  upon  volume  and  height. 

During  the  third  quarter  I  developed  code  for  sampling  uniformly  according  to  the  volume 
form  in  compact  subsets  of  two  dimensional  Minkowski  space,  and  in  cylinders  of  bounded 
height  with  flat  Lorentzian  metrics  (Minkowski  cylinders). 

During  the  fourth  quarter  I  used  PLEX  to  compute  the  persistent  homology  of  antichains 
in  posets  sampled  from  Minkowski  cylinders  and  showed  that  they  have  persistent  first 
Betti  number  over  a  long  range  of  scales. 

During  the  first  quarter  of  Year  2  of  the  grant  I  investigated  the  possibility  of  estimat¬ 
ing  another  persistent  topological  property,  dimension,  for  specific  metric  datasets  rele¬ 
vant  to  political  science.  I  described  preliminary  results  of  this  analysis  in  a  talk,  “Long 
range  dependencies  in  violent  events  timeseries",  in  the  Temporal  Quantitative  Analysis 
Panel  session  at  the  American  Political  Science  Association  Annual  Meeting ,  and  in  a  talk 
“Long-range  temporal  and  spatial  correlations  in  conflict”  in  the  session  on  Alternative 
Conceptions  of  State  Strength:  the  Role  of  State  Capacity  for  Development  and  Peace  at 
the  Sixth  Pan-European  International  Relations  Conference. 

During  the  second  quarter  I  applied  the  geometrical  methods  developed  for  Lorentzian  data 
to  probability  measures  on  finite  sample  spaces,  related  by  the  Bayesian  order  [6].  In  this 
case  the  homological  structure  is  delicate,  and  is  only  persistent  for  coarse  measurements 

[7]. 

During  the  third  quarter  we  began  investigating  the  possibility  of  analyzing  the  order 
complex  of  a  set  of  partially  ordered  data  within  the  framework  of  persistence  homology. 

During  the  fourth  quarter  we  concluded  that  for  geometrical  posets  there  was  little  useful 
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information  that  could  be  obtained  from  persistence  calculations  about  the  order  complex. 
We  also  began  investigating  the  possibility  of  analyzing  sets  of  symmetric  or  Hermitian 
matrices,  partially  ordered  by  positivity  or  majorization. 

During  the  first  quarter  of  Year  3  of  the  grant  we  developed  methods  for  simulating  Boolean 
test  data. 

During  the  second  quarter  we  demonstrated  that  simulated  Boolean  test  data  could  be 
generated  with  persistent  nonzero  first  and  second  Betti  numbers.  As  with  the  geometrical 
data,  the  calculations  were  implemented  with  PLEX. 

During  the  third  quarter  we  sampled  points  uniformly  from  a  singular  Minkowski  space 
(“pants”)  with  spatial  topology  changing  from  Sl  to  Sl  x  S1,  both  intrinsically,  and  also 
from  a  2  + 1  dimensional  Minkowski  space  in  which  the  pants  were  embedded.  We  showed 
that  persistent  homology  calculations  on  antichains  at  increasing  heights  in  the  poset 
capture  the  topology  change  [8]. 

During  the  fourth  quarter  we  began  formulating  our  partially  ordered  data  persistence 
constructions  in  a  category-theoretic  framework. 

During  the  first  quarter  of  Year  4  of  the  grant  we  continued  to  frame  our  results  category- 
theoretically. 

During  the  second  quarter  we  began  drafting  a  paper  on  our  results  from  this  project. 

During  the  third  and  final  quarter  we  continued  drafting  the  paper  describing  our  results 
from  this  project.  The  current  version  of  this  draft  is  attached  [9]. 

3.  Management  Summary 

3.1.  Personnel 

This  grant  supported  the  work  of  the  PI,  David  Meyer;  his  collaborator,  James  Pom- 
mersheim  (Reed  College);  and  his  graduate  students  at  UCSD,  Durdu  Giiney  and  Ben 
Wilson. 

4.  Publications  and  references 
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Introduction:  metric  versus  partially  ordered  data 
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Metric  and  non-metric  data 
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Thus  the  data  points  Ct  are  most  naturally  not  samples  from  a  (high 
dimensional)  metric  space,  but  are  rather  samples  from  the  partially 
ordered  set  of  subsets  of  T,  ordered  by  inclusion. 


Goal:  topological  analysis  of  partially  ordered  data 
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For  Q  C  P,  the  up-set  of  Q,  J+(Q )  =  {x  €  P  |  e  Q,  q  X  x}. 
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Open  covers  of  antichains 
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Nerves  of  open  covers 
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Unlike  the  metric  situation,  the  Rips  complex  is  also  more  convenient 
for  calculations  because  it  has  vertex  set  A  C  J+(A),  the  vertex  set 
for  the  nerve  complex. 


bO 

O  co 

O  — 
Q_  Q_ 

°  E 

+-»  L- 

03 

<u  x 

JZ  CD 


03 

+-» 

CO 

CL) 

"D 


03 

U 


a; 

E 

o 

a; 

bO 


-  v 

o  c 

+-»  c 

_  o 

CO  w 

X  a5 

£  12 

£Z  CO 

O  c 

u  o 

03  +J 

E  3 

M— 

C  CD 
—  CO 

<u  = 

CO  CO 

Z  ■- 
03  _-M 

10  -n 

<L)  CD 
co  +-> 
O  3 
Q.  Q. 

_  E 

wi  g 

O  bO 


-Z  <L> 

<  -O 


CO 

"D 


c 

03 


E 


CO 

a; 

Cl 

E 

03 

X 

CL) 


I- 


o 

co 

c 

a; 


o 

o 


E 

CO 

03 

CO 


< 


T 


■  • 


i 

CO 

O 

Cl 

o 

o' 

<u 

N 


CD 
> 
■4— > 

03 

bO 

CD 


co 


CD 

_X 

~a> 

u 

03 

CL 

CO 


< 


O 


CD 
> 
■4— > 

u 

CD 

CL 

CO 

CD 


CD 

> 


o 


Ll_ 


this  is  an  equivalence  relation  with  two  classes. 
One  is  designated  future-directed;  this  desig¬ 
nation  extends  to  null  vectors  in  the  closure 
of  the  set  of  future-directed  timelike  vectors. 


Posets  from  Lorentzian  manifolds 


to  = 

CO  3 
£Z 

o 


to 

u 

co 


cd 


CD 

T  .i 

■  ■  ■ 

CD  O 

>  -M 


U 


o 

cd 

co 


CD 
U 
C V 


u 

cd 

> 


~o 

cd 

+-» 

u 

CL) 


c 

CD  ■  = 

bo  "O 

ro  i_ 

+->  3 

>,  tf 


O  CD  '■*— 

O  5  ■o 

^  ,  c 
<C  ro 


H 


Yl  'I 

h  3s 

c" 

c/)  +-» 

O)  TO 
"O  -C 
CD  +J 
U  r- 

£  u 

Q.  3 

CO 

^  « 
ro  £ 

CO  ^ 

TO  u 

U  — 
TO 


UJ 

i_ 

O 


co 

n 

TO 

u  s* 

TO 

CO 


CD 

CD 


+3  "D 

-  fs 


C 

o 

i_ 

CD 

“O 


TO 


TO 

Q_ 

TO 

CO 

~co 

CD 

TD 

CD 

U 

CD 

Q_ 

_>, 

"to 

CO 

TO 

u 


B  ■ 
td  ~o 

<  g 


< 

r-H 

H 

"D 

< 

o 

•Q 

oj 


CD 

“D 

CD 

i_ 

TO 

CO 

CD 

CO 

CO 

CD 

u 

o 


CD 

"D  .E 

O 

Q_ 

e  -a 

V- 

O 


CD 


o 

> 

TO 

CO 

CD 

c 


CO 

TO 

U 

o 

■M 

CO 


U 


_  o 


CD 

TD 


U 

O  Q_ 

CO 


c 

CD 

■M 

U 


CD 

E 

CD 


CD 


I-  o 


c 

TO 

'n 

+-» 

c 

CD 

i_ 

O 

_J 

TO 

E 

o 

u- 

tuO 

c 

~Q_ 

E 

TO 

CO 

bb 

E 

c 

o 

"O 

c 

'  \_ 
CD 
TD 

!k_ 

TO 

o 

>, 

TO 

_Q 

CO 

— I 

TD 

CD 

TO 

U 

C 

TO 

-M 

CD 

_C 

_Q 

-M 

O 

TO 

c 

■M 

TO 

o 

TD 

bjO 

■M 

c 

CD 

CO 

'c 

O 

TO 

Q_ 

4— 1 

CD 

CD 

TD 

■ 

~o 

CO 

c 

,o 

O 

M— 

■ 

u 

C 

o 

CO 

TO 

E 
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Compute  persistent  homology  for  the  filtered  Rips  complex  R(Uj+(a))- 
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A  1-cycle  is  created,  and  then  persists  until  it  gets  filled  in  as  past 
lightcones  cover  a  little  more  than  1/3  of  the  initial  spacelike  S1. 
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But  there  are  also  natural  partial  orders  on  the  manifold  of  probability 
measures.  Majorization  is  one,  but  there  is  also  . . . 
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The  Bayesian  order  [Coecke  &  Martin  2002] 
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Partially  ordered  measure  data 
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Persistent  homology  for  Bayesian  measure  poset  data 
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Only  six  0-cycles  persist;  these  correspond  to  the  2-cells  of  the  barycen- 
tric  subdivision  of  the  2-simplex. 


Persistent  homology  for  coarse  Bayesian  measure  poset  data 
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Persistent  homology  for  coarse  Bayesian  measure  poset  data 
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Persistent  topology  depends  on  a  filtered  family  of  open  sets.  These 
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The  same  construction  gives  a  persistence  complex  for  any  subposet, 
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Introduction:  metric  versus  partially  ordered  data 
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Metric  and  non-metric  data 
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Thus  the  data  points  Rs  are  most  naturally  not  samples  from  a  (high 
dimensional)  metric  space,  but  are  rather  samples  from  the  partially 
ordered  set  of  subsets  of  T,  ordered  by  inclusion. 


Goal:  topological  analysis  of  partially  ordered  data 
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For  p  G  P,  the  height  of  p,  h(p),  is  the  length  of  a  longest  (maximal) 
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As  in  the  metric  case,  the  filtered  complexes  N{Uj+^a))  and  N(Up), 
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for  calculations  because  they  have  vertex  set  J^(A)  or  P ,  both  larger 
than  A,  the  vertex  set  for  the  Rips  complexes. 
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Finally,  suppose  that  when  a  student  understands  a  concept  s/he  is 
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A  1-cycle  is  created,  and  then  persists  until  it  gets  filled  in  as  past 
lightcones  cover  a  little  more  than  1/3  of  the  spacelike  S1. 
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lightcones  cover  a  little  more  than  1/3  of  the  spacelike  S1 . 


to 

•M 

c 

CO 

Q. 


c  ^ 
c  to 

<u  g 
O)  _§ 
CD  c 

eL2 

10 


CO 

03 


U 


o 

CL 


"O 

CD 


03 

O 

10 


u  (L) 

C  r- 


03  <D 

£  o 

“O  -M 
CD  ^ 
U  -Q 

“I  ^ 


CO 

Q. 

C 

CO 

‘n 

•M 

c 

<D 


Cl  - 


±  CD 


i 


CO  T— I 
03  + 


03 


CD 

_Q 


CN 


~o 

r<  CD 

“D 
-M  "O 
CD  CD 

_J  _Q 


O 

CO 

CD 
O 
03 
CL 

CO  CO 


03 

T 

CD 

+-> 

03 


bjO 


A  Lorentzian  pair  of  pants 
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Persistent  topology  depends  on  a  filtered  family  of  open  sets.  These 
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Abstract 


In  this  paper,  we  study  the  topological  properties  of  posets  using  the  tool  of  persistent 
homology.  We  introduce  a  construction,  given  a  poset,  which  results  in  a  filtered  simp  lie  ial 
complex.  We  investigate  this  construction  on  posets  sampled  from  Lorentzian  manifolds. 


1  Introduction 

2  Preliminaries 

We  now  remind  the  reader  of  some  elementary  facts  about  posets.  Let  P  be  a  finite  set  throughout. 

Definition  2.1.  A  partially  ordered  set  or  poset  (P,  is  a  set  P  with  a  partial  order  ■<.  If  two 
elements  a,  b  E  P  are  not  related,  we  write  a  /  6. 

A  partial  order  is  a  binary  relation  satisfying  reflexivity,  antisymmetry,  and  transitivity.  Through 
out  this  paper  we  make  use  of  special  subsets  of  posets  in  which  no  two  elements  are  related. 

Definition  2.2.  An  antichain  A  is  a  subset  of  a  poset  (P,  such  that  the  following  holds: 

•  For  any  two  a,  b  E  A,  a  /  b. 

Example:  Let  n  €  N.  Let  Sn  denote  the  following  poset.  The  elements  of  Sn  are  the  ax  and  bt 
for  1  <  i  <  n  with  the  relations  a*  <  bj  for  i  ^  j. 


•a  i 


•a2 


•U3 


•a4 


The  sets  of  {ai, . . . ,  an}  and  {6i, . . . ,  6n}  are  two  antichains  in  Sn. 

When  discussing  posets,  it  is  often  easiest  to  describe  immediate  relations  and  allow  all  other 
relations  to  be  implied  by  transitivity. 


Definition  2.3.  Given  a  finite  binary  relation  (P,  ^b),  the  transitive  closure  (P,  ^p)  is  the  minimal 
transitive  order  which  extends  ; <b •  That  is  for  a,b  €  B,  a  ~<p  b  if  there  exists  ci, . . .  ,c*  6  P  such 
that  a  Cl  ~<B  C2  diB  •  *  •  diB  ck  ^±B  b. 


The  transitive  closure  of  a  partially  ordered  set  is  the  partially  ordered  set  itself.  Given  an 
antichain  A  in  a  poset  P,  we  can  define  a  height  function  which  measures  the  distance  to  the  given 
antichain.  We  will  denote  a  path  from  p  to  q  by  elements  po,  •  •  •  ,Pn  £  P  where  p  =  po  ^  p\  ■<  •  •  •  ■< 
pn  =  q.  The  length  of  the  path,  £{po,  •  •  •  ,Pn)^  is  n.  Let 


d(p,  q) 


min 


paths  (p=poi-,Pn=q)€TP 


t{PO,---,Pn) 


If  there  are  no  paths  from  p  to  g,  this  distance  is  oo.  Finally  we  define  our  distance  function  as 


d(p,q)  =  min  (d(p,q),d(q,p)) 


Now  we  can  define  our  height  function.  Let  A  be  an  inextendable  antichain.  Define  :  P  — +  Z  as 
follows: 


hj\(p)  =  min{d(a,p)  :  a  G  A} 
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This  function  is  well  defined  because  A  is  inextendable.  Suppose  hA(p)  =  oo  for  some  p  G  P. 
Then  there  are  no  paths  starting  in  A  and  ending  at  p,  nor  paths  starting  at  p  and  ending  at  A. 
Thus  p  is  incomparable  with  all  elements  of  A.  So  A  can  be  extended  by  p.  Two  noted  inextendable 
antichains  are  the  sets  Am  and  Am  of  minimal  and  maximal  elements  of  the  poset  respectively.  Let 
hm  and  Hm  denote  the  corresponding  height  functions. 

Proposition  2.4.  Aj  =  {p  G  P\hm(p)  =  j}  U  {p  G  P\hm(p)  <  j  and  p  is  maximal  }  is  an  inex¬ 
tendable  antichain  in  the  poset  for  any  j  G  Z. 

Example: 


#  •  • 

IX  XI 


ft  ft  ft 


2 

1 

0 


•  1 

(P<)  hA(p) 

3  Construction 

We  now  discuss  the  construction  of  a  simplicial  complex  given  a  poset  and  antichain  pair.  Let 
FinSimp  denote  the  category  of  simplicial  complexes  and  simplicial  maps.  Let  PA  denote  the 
category  of  pairs  (P,  A)  of  finite  posets  P  and  antichains  A  C  P.  For  objects  (P,  A),(R,  B)  E  PA , 
Hom((P,  A),(R,  B))  is  the  set  of  order  preserving  maps  f  :  P  —>  R  with  f(A)  C  B.  Given  an 
object  (P,  A)  of  PA ,  with  A  =  {ai, . . . ,  aw},  one  can  construct  an  abstract  simplicial  complex.  The 
simplicial  complex  will  have  vertices  a\, . . .  ,ax-  What  we  need  is  a  rule  dictating  whether  or  not 
to  assigning  a  fc-simplex  to  fc  + 1  points.  This  assignment  will  be  determined  by  other  points  in  the 
poset. 

Definition  3.1.  Let  p  G  P.  Define  I(p)  =  {q  G  P\p  ^  q  or  q  X  p).  Define  Av  =  I(p)  Pi  A. 

The  set  I(p)  is  the  set  of  elements  related  to  p  in  the  poset.  Each  Ap  is  the  set  of  elements  in  A 
related  to  p,  which  gives  us  information  about  the  relative  closeness  of  elements  in  the  antichain. 
That  is,  two  elements  of  A  which  appear  together  in  some  Ap  are  closer  together  than  two  elements 
of  A  which  do  not.  We  use  this  information  to  construct  a  functor  A  :PA  — *  FinSimp.  Given  a  poset 
antichain  pair  (P,  A),  the  vertices  of  the  simplicial  complex  A(P,  A)  are  {ai, . . .  the  elements 

of  the  antichain.  The  vertices  a^, . . . ,  aik  form  a  A;-simplex  if  there  is  a  p  G  P  with  alQ, . . . ,  X  p. 
Alternatively,  if  Ap  =  {a;0, . . . ,  a*fc},  we  add  the  fc-simplex  [af0, . . . ,  aij  to  the  abstract  simplicial 
complex  A(P,  A).  This  creates  a  well  defined  simplicial  complex  as  the  condition  will  include  all 
subsimplices  of  a  given  simplex  by  definition. 

Example:  For  the  poset  Sn  and  antichain  A  =  {ai, . . .  ,an},  A (Sn,A)  is  the  boundary  of  the 
standard  n  —  1  simplex.  This  simplicial  complex  is  homotopy  equivalent  to  Sn'2. 

Lemma  3.2.  This  construction  only  depends  on  the  minimal  and  maximal  elements  of  P. 

Proof  If  p  G  P  is  neither  minimal  nor  maximal,  then  there  is  a  q  e  P  with  p  <  q  or  q  <  p.  Assume 
p  <  q.  Then  Ap  C  Aq.  Thus  the  simplex  associated  to  Ap  is  a  face  of  Aq.  □ 
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If  P  has  a  unique  minimal  or  maximal  element,  the  resulting  complex  A(P,  A)  is  a  |A|-simplex 
for  any  antichain  A  in  P. 

It  is  enough  to  define  a  map  /  :  K  — ►  L  of  simplicial  complexes  by  giving  a  map  of  vertices  so 
long  as  the  following  holds.  If  [tq, . . . ,  v*]  is  a  simplex  in  K ,  then  [f(v i), . . . ,  f{vk)]  must  be  a  simplex 
in  L.  We  now  show  this  holds  for  our  constructed  complex.  Let  /  :  (P,  A)  — *  (P,  B)  be  a  morphism 
in  PA ,  that  is  an  order  preserving  map  with  f(A)  C  B.  Define  A (/)  :  A (PA)  — ►  A {R,B)  in  the 
following  way.  Let  A  =  {ai, . . .  ,  a;v}.  Set  A(/)(aj)  =  /(at).  This  defines  a  map  on  the  vertices  of 
A(P,  A).  We  need  only  check  that  simplices  are  mapped  to  simplices.  Suppose  [an, . . . ,  a{k ]  is  a 
simplex  in  A(P,  A).  Then  there  is  a  p  G  P  such  that  {aq, . . .  C  Ap.  Then  /(a*)  G  Bf for 
all  z,  by  order  preservation.  So  the  simplex  [/(an),.. .  ,f{a>ik )]  is  in  A (P,P)  and  the  map  is  well 
defined. 

There  is  a  functor  T  :FinSimp^PA  going  the  other  direction,  similar  to  A.  Let  K  be  a  finite 
simplicial  complex.  F(K)  is  a  poset  P  given  by  the  following  construction:  Let  Ki  denote  the  set 
of  z-dimensional  simplices  of  K.  The  constructed  poset  will  have  ^2\K{\  elements.  By  inducting  on 
dimension,  we  can  create  a  poset  by  defining  immediate  relations  and  taking  the  transitive  closure. 
For  i  =  0,  let  Kq  =  {r>o>* . .  ,  v/v}  denote  the  set  of  vertices  of  K.  Add  an  element  to  P  for  each 
element  of  K§.  For  i  =  1,  add  an  element  to  P  for  each  element  of  K\.  Then,  for  each  element 
e*.  =  °f  add  relations  ^  and  ^  £k-  That  is,  all  of  the  subsimplices  of  e ^ 

are  less  than  e*.  We  then  proceed  inductively.  For  each  i  €  N,  we  add  an  element  to  P  for  each 
element  in  Ki ,  with  relations  to  all  of  its  subsimplices.  We  set  A  to  be  the  set  of  minimal  elements 
of  the  constructed  P. 

Example:  The  circle  with  three  vertices  and  three  edges  is  mapped  to  the  crown  poset  with  six 
elements  of  height  2  and  antichain  the  set  of  minimal  elements  under  T. 


Proposition  3.3.  A  o  T  :FinSimp  — >  FinSimp  =  Idf  in$imp 

4  Persistence 

The  following  definition  of  persistence  comes  from  [Carlsson08] 

Definition  4.1.  Let  C  be  a  category  and  P  be  a  poset.  We  can  regard  P  as  a  category  V  in  the 
usual  way,  z.e.,  with  object  set  P  and  a  unique  morphism  from  x  to  y  whenever  x  <y.  Then  by  a 
P-Persistent  Object  we  mean  a  functor  :  V  — ►  C. 

The  paper  by  Carlsson  gives  examples  of  N  and  R  persistent  objects.  A  natural  category  of 
study  is  the  category  of  topological  spaces.  Morse  theory  describes  a  continuous  filtration  of  a 
topological  space  given  by  a  Morse  function.  This  filtration  of  the  topological  space  produces  an 
R-persistent  object  in  the  category  of  topological  spaces. 

Definition  4.2.  An  N -filtration  of  a  space  A  is  a  sequence  X()  C  A1  C  A2  C  . . .  of  subspaces  of 
A  whose  union  is  A. 

Given  a  poset  P  and  antichain  A,  wre  will  now  describe  an  N-filtration  on  the  associated  simplicial 
complex  A(P,  A).  Suppose  P°  C  •  •  •  C  PA/  =  P  is  an  increasing  sequence  of  sub-posets  of  P.  Then 
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we  can  define  a  filtration  on  A (P,  A)  by  A (P,  A)i  =  A (P-7,  >1  fi  P-7)  for  each  j  £  N.  That  is  to  say, 
the  jth  step  in  the  filtration  is  the  simplicial  complex  obtained  from  the  construction  associated  to 
the  poset  P-7  with  antichain  P-7  C\A.  Given  a  height  function  ft  as  described  above,  we  get  a  natural 
sequence  of  sub-posets.  We  can  set  P7  =  {p  £  P  :  ft(p)  <  j},  filtering  the  poset  by  its  distance 
from  the  antichain. 

Theorem  4.3  (Realization).  Let  K  be  a  finite  simplicial  complex.  Let  {F1}  be  an  N -filtration  of 
the  simplicial  complex  with  K(0 )  =  F°.  There  is  a  poset  P  and  an  antichain  A  C  P  such  that 
A(P,A)  =  K .  Moreover,  the  shortest  path  height  function  h  provides  the  same  filtration  as  given 
by  {F*'}. 

Proof.  Suppose  K =  {vo, . . . ,  vn}  and  let  P  be  an  empty  poset.  Let  Me  N  be  such  that 
pi  _  pM  £or  ajj  j  >  jrf  por  each  simpiex  ailmm,ik  =  [vq , . . . ,  Vik)  in  F1,  add  an  element  a}x  ifc  to 
P.  Then  add  the  relations  v®.  •<  af  for  j  =  1 . . .  k.  For  each  simplex  <Jilm..ik  —  ,  v<fc] 

in  F2,  add  an  element  of  to  P.  Then  add  the  relations  ah  -<  a 2  4,  for  j  =  l . .  .k.  Proceed 

*i  •••*«  *i*«-*fc 

inductively.  For  a  simplex  =  [vq, . . .  ,Vik]  in  FL  add  an  element  o\x  ik  to  P.  Then  add 

the  relations  o\~l  <  o\x  ifc  for  j  =  1 . . .  k.  Finally,  take  P  to  be  the  transitive  closure  of  P.  This 

construction  will  result  in  chains  v®  ■<  o }  •<  •  •  •  ■<  of!~l  <  erf'1  for  each  vertex  present  in  K ,  with 
additional  elements  and  relations  in  P  for  each  simplex  present  in  K 


-M 

7o 


* 


•_A/  —  1  •  M  1 

(X0  <T1 


°JV-1 


.M- 

N- 


•_M  —  1 


2 

iV-1 


2 

N 


•„! 


1 

N-l 


1 

jV 


•„o 


Now  to  show  this  construction  has  the  desired  properties.  Let  A  =  {p  £  P  :fiq  £  P  with  <?  <  g}, 
the  set  of  minimal  elements  of  P.  Let  ft  :  P  — *  N  be  the  length  of  the  shortest  path  to  A  function. 
Let  P-7  =  {p  £  P  :  ft(p)  <  j}  denote  the  filtration  of  P  given  by  this  height  function.  Bv 
construction,  A(P°,  P°)  =  K ^  =  F°.  Now  for  a  general  j  £  N  with  j  <  M.  Let  a  =  . . .  ,Vik] 

be  a  simplex  in  F7.  Then  -fc  £  P.  Moreover,  ft(^A  lk)  <  j  because  the  chain  of  relations 
v®  ^  o\  <  •  *  •  ^  ^  ifc  exists  for  each  t  =  1 ...  ft.  So  the  shortest  path  to  A  has  length  at 

most  j  and  o\x  ik  £  P*7.  Thus  o  is  a  simplex  in  A (P*7,  P^H^l)  and  we  have  that  F-7  C  A(P7,  P^D^l). 
Now  let  rr  be  a  simplex  in  A (P-7,  P-7  D  ^4).  The  constructed  simplicial  complex  A (P7,  P-7  fl  ^4)  only 
depends  on  the  maximal  elements  of  PJ .  So  is  constructed  due  to  an  element  lk  £  P-7.  This 
element  was  added  in  correspondence  with  the  simplex  [v^, . . . ,  Vifc]  £  F*7.  By  construction,  the 
element  o\x  ik  £  PJ  only  adds  [v^ , . . . ,  Vik]  and  its  sub-simplices  to  A (P-7,  P7  fl  >1).  So  cr,  must  be 
[t?n, . . .  ,V{k]  or  one  of  its  sub-simplices.  Thus  A (P-7,  P-7  fl  ^4)  =  F*7.  □ 


Example:  Given  the  2-simplex  filtered  in  the  following  way,  the  above  construction  produces 
the  following  poset. 
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This  construction  is  not  unique.  There  are  many  other  posets  that  will  create  the  same  filtered 
simplicial  complex. 

5  Category 

Let  P  =  {pi, . . . . p\ }  be  a  finite  poset  and  A  an  antichain.  Let  By  denote  the  boolean  algebra 
on  N  elements.  The  category,  Bn,  has  a  partial  order  given  by  inclusion.  There  is  a  functor  from 
Bn  to  VA  which  takes  an  element  {ii  > . . . ,  u }  of  Bn  to  the  subposet  {pri , . . . , pik }  and  which  takes 
morphisms  in  Bn  to  inclusions  of  posets.  The  category  Bn  is  a  directed  set. 

Proposition  5.1.  lim  I  =  {1, . . . ,  N) 

i7bn 

Proof.  {1, . . . ,  AT}  is  cofinal  in  Bn  because  {1, . . . ,  N}  is  the  unique  maximal  element  of  Bn •  □ 

For  a  fixed  antichain,  we  can  define  a  category  Cp(A).  The  objects  of  this  category  are  A (5,  S  fl 
^4)  for  each  S  C  P,  denoted  A 5.  That  is,  this  category  is  the  image  of  the  above  functor.  The 
morphisms  are  defined  as 


HomcP(i4)(A5l,As2) 


l  :  ASl  -  AS2  if  Si  C  S2, 
0  otherwise. 


This  category  is  a  directed  set.  The  object  A p  6  Obj(Cp(A))  is  cofinal. 


Proposition  5.2.  limA^  =  Ap 


With  definition  4.1  in  mind,  define  :  Bn  —*  Cp(A)  to  be  ^(Pi)  =  A px.  Thus  the  whole 
category  Cp{A)  is  a  Bn- persistent  object. 

Proposition  5.3.  Let  P  be  a  poset  and  :  V  — ►  C  be  a  P-persistent  object .  That  is ,  we  have  a 
family  of  objects  {cx}*€P  C  Obj(C)  together  with  morphisms  <j)Xy  '  cx  —+  Cy  whenever  x  ■<  y  such 
that  <fyz  o  (pxy  —  (f>xz  whenever  x  <y  <  z.  Suppose  po  ^  P\  r<  •  •  •  d:  PN  is  a  chain  in  P.  Then  there 
is  an  N -persistent  object  in  C  given  by  the  functor 
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$ 


POPl-PN 


(i)  = 


if  0  <  i  <  N 
if  i  >  N 


Proposition  5.4.  If  T  :  C  — >  V  is  a  functor  and  $  :  V  —*  C  is  a  P-persistent  object ,  then 
T  o  4>  :  V  — i  V  is  a  P-persistent  object. 


6  Simplicial  Homology 


Once  we  have  a  simplicial  complex  built,  a  natural  set  of  invariant  to  compute  are  the  Betti  numbers. 
We  give  a  brief  description  of  the  Betti  numbers  here.  For  a  more  detailed  reference,  see  [Munkres]. 

Definition  6.1.  The  nth  chain  group  of  an  abstract  simplicial  complex  A',  denoted  Cn(K),  is  the 
free  abelian  group  generated  by  the  oriented  n-simplices  of  A.  That  is  [a]  =  [r]  if  a  =  r  and  a  and 
r  are  oppositely  oriented. 

Definition  6.2.  The  simplicial  chain  complex  associated  to  K  is  the  chain  complex  •  •  •  — ►  C2(K)  — ^ 
C\{K)  — +  Co(A')  ^  0  — ►  . . .  with  boundary  operator 

n 

d„(K,,...,Vin])  =  £(-!)>*,. 

j=o 

where  [• Vi0 , . . .  . . .  ,Vjn]  denotes  the  n  —  1  simplex  formed  by  omitting 

It  is  a  standard  exercise  to  show  that  dx  o  dx+\  =  0  for  each  i  £  N.  This  makes  it  possible  to 
define  our  homology  groups. 

Definition  6.3.  The  nth  homology  group ,  Hn(K),  is  given  by  Hn(K)  =  keidn/  im  dn-fl- 
Definition  6.4.  The  nth  Betti  number ,  f3n(K)y  is  defined  as  i3n(K)  =  rankz (Free (Hn {K))) 
Example:  The  circle 


( 


Hn(Sl) 


if  n  =  0, 1 
otherwise 


Hn  (51) 


1  if  n  =  0, 1 
0  otherwise 


7  Dimensionality 

Let  (P,  A)  be  a  finite  poset  antichain  pair.  Let  A(P,  A)  be  the  associated  simplicial  complex. 

Definition  7.1.  The  dimension  of  a  poset  P  is  the  minimum  number  of  total  orders  whose  inter¬ 
sections  give  the  partial  order  on  P 

Example:  The  poset  Sn  has  dimension  n. 

Proposition  7.2.  If  Q  C  P  is  an  induced  subposet ,  then  dim  Q  <  dimP. 
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Let  A"  be  a  finite  simplicial  complex.  Let  TV  denote  the  number  of  vertices  present  in  X.  Then 
A"  is  isomorphic  to  a  subcomplex  of  A  v,  the  standard  JV-simplex.  So  the  polytope  of  X ,  |Ar|, 
embeds  in  R'v+1. 

Definition  7.3.  Let  a  be  a  simplex  in  A\  Let  V  =  . . . ,  e{k}  be  the  vertices  of  a.  Let 

7 Ta  :  RA+1  — >  R*v+1  be  the  projection  map  onto  the  coordinates  of  a.  Let  7r£  =  Id  —  7 ra.  That  is 
projection  onto  the  coordinates  not  in  a. 

If  r  E  R+,  Define  U{a,  r)  =  {x  :  x  E  R%and  £  7ra(x)t  >  1  —  r}. 

Define  V(a,r)  =  U(a,r)  D  \X\. 

Va%r  is  open  in  X  since  the  topology  on  A”  is  the  same  as  the  induced  topology  from  R'V+1,  and 
r)  is  open  in  R^+1. 

Note:  If  r  <  1,  V(a,r)  can’t  contain  points  in  a  simplex  which  does  not  intersect  a.  If  r  < 
then  V (er,  r)  D  V(r,  r)  =  0  whenever  a  D  r  =  0. 

Lemma  7.4.  If  a  E  X  and  r  <  1,  V(a,  r)  is  homotopy  equivalent  to  a. 

Proof.  Define  H(x,  t)  :  V(a,r)  — *  V(a,i')  by 

H(x,  t)  =  (1  -  t)x  +  t-=£z^rr 

/  >  'Ka(X)i 

Note  this  is  well  defined  since  7r <j(x)i  >  1  —  r  ^  0. 


7Tcr(ar)i 


',n<r(x)i 

-  (!  -  0  Y  xi  +  73  Y  **(*)* 

^  Tra(x)i  *—< 

=  (1  —  t)  +  t 

=  1 


So  the  homotopy  remains  in  A  % .  The  sum  of  the  coordinates  is  always  1,  and  if  x,  =  0  for  any 
z,  then  H(x,  t)i  =  0  for  all  t.  Thus,  if  x  6  r  for  some  r  E  A",  then  H(x ,  t)  €  r  for  all  £  €  [0, 1].  The 
homotopy  fixes  points  in  a.  So  a  is  a  deformation  retract  of  V(cr ,  r).  □ 

Since  a  is  contractible  in  |AT|,  V(cr,  r)  is  contractible  in  \X\. 

Lemma  7.5.  If  <r,  r  E  X,  then  V(ar ,  n)  D  V*(t,  r2)  C  V(a  Dr,  r\  +  r2). 

Proof.  Let  x  E  V(cry  ri)  n  V{r,  r2).  Then 

y^7rff(x)j  >  l-n, £>(*)<  >  1_r2 

and 

Since  the  coordinates  sum  to  1,  we  also  have  that 


Y  <  rl  ’  Y  nr(x)i  <  r2 


So 


<  r2 
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Y  n°nr(x)  =  Y  "SC*)  +  Y  *?(*)  ~  Y  "SurW 

X  *5nr(*)  <  ft  +  r2  ~  Y  *5ur(*) 

53*5nr(*)  <  n  +  r2 


□ 

If  ri,r2  <  1  and  <j  D  r  =  0,  then  V(cr,ri)  fl  V(r,  7-2)  =  0.  If  n  +  r2  <  then  V(cr,  ri)  n 

V^r,  r2)  =  {x  G  \X\  :  53  >  1  —  r\  and  53  frr(x)i  >  1  —  r 2}  is  contractible  in  \X\.  The  above 

homotopy  gives  the  retraction  onto  cr  D  r,  which  is  contractible.  In  fact,  if  a\  fl  *  •  •  fl  an  ^  0  and 
ri  +  T2  4 - h  rn  <  1,  then  F((7i,ri)  fl  ^(<72,1*2)  fl  •  •  •  fl  V(crn,rn)  is  contractible. 

Lemma  7.6.  If  X  is  a  simplicial  complex  with  at  most  K  facets,  then  Hn(X)  =  0  for  n>  K  —  2. 

Proof  Embed  \X\  in  Aa  C  Ka+1  as  above.  Let  M  =  {cr  G  X  :  cr  <£  r  for  r  ^  a  G  X},  the  facets 
of  X.  By  assumption  |A/|  <  K .  Let  (9  =  {I/(cr,  :  cr  G  A/}.  Then  0  is  an  open  cover  of  A" 

since  every  x  G  X  is  in  some  facet  of  X.  For  any  cr  G  A/,  V(cr ,  ^)  is  contractible.  By  the  above 

lemma,  and  the  fact  that  the  intersection  of  two  simplices  is  a  face  of  both  simplices,  we  see  that 

the  intersection  of  two  such  open  sets,  if  nonempty,  is  also  contractible.  For  crj, . . .  ,cr *  €  M, 

n  V(<T2, 1) n •••  n  vfo,  A)  c  F(r, 

where  r  =  crj  fl  •  •  •  fl  cr^.  Since  A:  <  A",  ^  <  1  and  V(r,  is  contractible  in  X.  Since  O  is  a  good 
cover  of  |A|,  the  nerve  of  0,  N(O)  has  the  same  homology  as  |A|.  Since  this  cover  has  only  K 
elements,  all  K  +  1  intersections  are  empty  and  there  is  only  one  K  intersection  to  consider.  If 
all  of  the  facets  of  X  intersect  in  a  common  simplex,  then  the  K- wise  intersection  is  nonempty. 
Otherwise,  the  nonempty  K  wise  intersection  results  in  a  single  K  —  1  dimensional  simplex  in  the 
nerve.  The  boundary  of  this  simplicial  complex  is  nonempty.  In  either  case,  the  K  —  2th  homology 
of  the  nerve  simplicial  complex  is  trivial.  For  dimension  reasons,  all  higher  homology  groups  are 
trivial  as  well.  □ 

Theorem  7.7.  Suppose  Hn(A(P,A))  ^  0.  Then  dim  P  >  [(n  +  2)/2].  If  A  is  a  minimal  or 
maximal  antichain ,  then  dim  P  >  (n  +  2) 

Proof  By  the  above  lemma,  A (P,  A)  must  have  at  least  n  +  2  facets.  Let  V  =  {i>i, . . .  be  the 
vertices  present  in  A(P,  ^4)  and  /1, . . . ,  ft  be  the  facets  present  in  S.  By  construction,  each  of  these 
facets  is  present  due  to  an  element  p  G  P.  Let  I(V)  be  the  past  and  future  of  V.  To  show  at  least 
n  +  2  distinct  elements  of  M  are  needed  to  create  S.  Suppose  M  =  generates  all 

of  the  faces  of  S.  Then  S  can  be  thought  of  as  a  nonzero  element  of  the  nth  homology  group  of 
A(Af  U  A,  A).  The  boundary  of  S  is  still  empty.  If  5  was  not  the  boundary  of  an  n  +  1  chain  in 
A(P,  A))  then  it  will  not  be  A(A/  U  ^4,^4)  either,  since  there  are  fewer  such  chains.  By  the  above 
lemma,  K  >  n  +  2. 

We  now  list  total  orders  on  {m*}  U{  nt}  which  must  be  present  in  any  realizer  of  this  partial 
order. 

This  subposet  consists  of  v\, . . . ,  1,7+2  which  are  unrelated  and  mi, ... ,  mn+2.  For  each  2,  is 
related  a  unique  set  of  vertices,  and  either  preceeds  all  of  them,  or  follows  all  of  them.  Split  the 
mt ’s  into  two  groups.  One  in  which  the  ml1s  preceed  vertices,  call  it  A/iower  and  one  where  the  m^’s 
follow  vertices  A/upper. 
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For  each  mj  G  Miower  there  must  be  a  total  order  <;  on  V  with  m*  <;  vj  for  all  Vj  G  Ami  n  V 
and  Vj  <i  mi  for  all  Vj  Ami  Pi  V.  Since  mi  has  a  distinct  set  Ami  of  vertices,  we  must  have  a  total 
order  for  each  mi  G  A/iower- 

Similarly,  we  need  a  total  ordering  wrhich  gives  the  relations  for  each  mi  G  A/Upper-  These  upper 
and  lower  mf  s  can  be  related,  so  the  total  orderings  of  one  may  give  the  total  orderings  of  the 
other.  So  we  need  at  least  max  |Miower|,  |MupPer|  total  orderings  in  any  realizer  to  induce  the  proper 
relations.  By  the  pigeon  hole  principle,  this  max  is  at  least  [|M|/2]  =  [(n  +  2)/2].  Thus  we  get 
that  dimP  >  \(n  +  2)/2].  If  A  is  a  minimal  or  maximal  antichain,  then  the  above  max  is  |M|  and 
we  have  that  dimP  >  (n  +  2). 

□ 

This  bound  is  tight.  For  the  poset  Sn+2  and  antichain  A  =  {ai,.  . .,an+2}>  A(Sn+2yA)  ~  Sn. 
This  poset  has  dimension  n  +  2  and  /?n(5n)  =  1. 


8  Persistent  Homology 


We  have  described,  for  a  given  poset  P,  antichain  A  and  height  function  h,  a  filtered  simplicial 
complex  A(P,  A).  One  can  compute  the  homology  of  the  simplicial  complex  A(P,  A)  at  any  given 
step  of  the  filtration  to  obtain  some  information  about  the  original  poset.  A  more  useful  approach 
is  to  compute  the  homology  at  every  step  of  the  filtration  as  a  summary  of  information  about  the 
original  poset. 

Let  T  be  a  totally  ordered  set.  Let  R  be  a  ring.  Now  to  describe  a  special  type  of  persistent 
object,  one  whose  associated  poset  is  totally  ordered  and  target  category  is  the  category  Rmod. 

Definition  8.1.  A  persistence  module  parameterized  by  T  over  R  is  a  family  of  P-modules  {Mt}teT 
with  P-module  homomorphism  {<£*,*' }(t,t')€TxT  su°h  that  iptj  o  —  {Pt>tn  whenever  t  <t!  <  tn . 

When  working  over  a  field  P,  finite  type  persistence  modules  parameterized  by  Z  have  a  nice 
decomposition,  as  formulated  in  [CZ04].  We  now  describe  this  decomposition.  Let  I  =  (m,n)  c 
N  U  {oo}  be  an  interval.  Define  the  persistence  module  Q{I)  parameterized  by  N  over  F  to  be 


Q(i)t  = 


if  tel, 
otherwise. 


Definition  8.2,  A  persistence  module  parameterized  by  N  over  F,  { Mn}n is  said  to  be  of  finite 
type  if 


1.  Mn  is  a  finite  dimensional  F  vector  space  for  each  n  G  N 


2.  There  exists  an  N  such  that  for  all  t  >  N,  is  an  isomorphism. 

Proposit  ion  8.3.  A  persistence  module  parameterized  by  N  over  F  of  finite  type  is  isomorphic  to 
one  of  the  following  type 

N 

0  Q(h) 

5=1 

where  Is  C  N  U  {oo}  is  an  interval  for  each  s  G  N. 

In  short,  to  encode  the  information  to  describe  such  a  persistence  module,  we  only  need  to  list 
the  relevant  intervals.  Hence  the  definition 

Definition  8.4.  A  barcode  is  a  collection  of  intervals  {IS}^=1,  where  each  interval  Is  is  bounded 
below. 
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By  the  definition  of  filtration,  we  have  an  inclusion  map  i :  A (P,Ay  *— ►  A(P,  A)k  for  all  j  <  k. 
So  we  have  a  N-persistent  object  in  the  category  of  simplicial  complex  whenever  we  have  an  N 
filtration.  By  the  functoriality  of  Hn(— ),  wTe  then  get  a  persistence  module  parameterized  by  N 
over  Z. 

Example:  The  circle  filtered  in  the  following  way  produces  the  following  persistence  modules 
and  barcodes. 


H0  =  Z3  Z  Z 

Hi  =0  0  Z 


The  circle  filter  in  this  manner  has  Betti  0  barcode 

and  Betti  1  barcode  INSERT  BARCODE  PICTURE  HERE 

Example:  Let  P  be  the  crown  poset  of  height  2  on  3  elements.  Let  A  be  the  antichain  of 
minimal  elements  of  P.  Let  h  be  the  height  function  such  that  the  minimal  elements  have  height 
0  and  the  maximal  elements  have  height  1. 


The  poset  filtered  by  height  is  then 


•  height  1 

•  height  0 


•  •  • 


•  •  •  •  •  • 

This  filtration  on  the  poset  results  in  the  filtered  simplicial  complex 

•  • 

/ 

\ 

•  « 

9  Sampling  from  a  Lorentzian  Manifold 

One  can  take  a  finite  sample  of  a  Lorentzian  manifold  and  create  a  poset  using  the  causal  structure 
of  the  manifold.  This  poset  should  contain  some  of  the  topological  information  about  the  manifold. 
By  picking  antichains  in  the  created  poset  (analogous  to  picking  space-like  hypersurfaces  of  the 
original  manifold)  we  can  investigate  topological  properties  of  the  poset. 
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10  Computed  Examples 


Let  (M,  g)  be  a  Lorentzian  manifold.  Let  X  C  M  be  a  finite  collection  of  points.  The  set  A”  can  be 
given  a  partial  order  using  the  causal  structure  of  the  manifold.  For  two  points  x,y  €  X, 

we  say  that  x  ^(Myg)  V  if  there  is  a  future-directed  time- like  curve  from  x  to  y  in  M. 

Cylinder  In  this  example,  we  attempt  to  recover  the  space-like  topology  of  the  Lorentzian 
cylinder.  The  cylinder  is  described  as  an  identification  space  of  the  unit  square,  and  is  given  the 
flat  metric.  1000  points  were  sampled  the  unit  square.  The  set  was  given  the  partial  order  based  on 
the  (1+1)  Lorentzian  metric  on  the  cylinder.  The  elements  were  assigned  a  height  based  on  their 
shortest  distance  to  the  minimal  elements  of  the  poset.  antichains  in  the  poset  were  chosen  to  be 
sets  of  elements  at  the  same  height. 

Example:  Height  0,  the  minimal  elements  of  the  poset.  The  filtered  construction  resulted  in 
the  following  barcodes 

Pair  of  Pants  In  this  example,  we  attempt  to  recover  the  space-like  topology  of  the  Lorentzian 
pair  pants.  The  pair  of  pants  were  created  using  an  identification  space  of  the  unit  square. 


The  pair  of  pants  minus  the  singularity  was  given  the  flat  metric.  2000  points  were  sampled 
from  the  space  and  given  the  partial  order  based  on  the  Lorentzian  metric.  Elements  of  the  poset 
were  assigned  a  height  based  on  their  shortest  distance  to  the  minimal  elements  of  the  poset. 
antichains  in  the  poset  were  chosen  as  sets  of  elements  of  the  same  height.  The  resulting  barcodes 
were  dependent  on  the  height  of  the  antichain,  antichains  taken  below  the  singularity  resulted  in 
single  persistent  first  homology  elements,  whereas  antichains  taken  above  the  singularity  resulted 
in  a  pair  of  persistent  first  homology  elements. 

INSERT  BARCODE  PICTURES  HERE 
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